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This meeting produces another evidence that present parallel computers are (a) real 
instruments of computational physics, (b) largely in the hands of still-pioneers, (c) effi- 
ciently promoted by basic research groups with large-scale computationalneeds. Progress 
in parallel computing is carried by two types of such groups, that either follow the build- 
it-yourself or the early-use strategies. In this contribution, we describe, as an example 
to the second approach, the Wuppertal university pilot project in applied parallel com- 
puting. We report in particular about one of our key applications in theoretical particle 
physics on the Connection Machine CM-2: a high statistics computer experiment to de- 
termine the static quark-antiquark potential from quenched quantum chromodynamics. 
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1. Outline 

In section 2 we present a brief expose on quantum chromodynamics to non-experts, 
in order to illuminate our motivation to invest for considerable efforts in large-scale 
computing. We will also shortly mention the Wuppertal activities on parallel com- 
puting, based on the pilot installations of the Connection Machines CM-2 and CM-5. 
In section 3, we will illustrate that parallel computers have matured to real research 
tools in numerical quantum field theory. We will examplify this by presenting the 
state-of-the-art of Creutz' classical 1799 computer experiment, employing present 
hardware and software techniques in parallel computing. Section 4 contains a brief 
physics discussion of our results. 

*Talk presented by K. Schilling at the "Workshop on Large Scale Computational Physics on 
Massively Parallel Computers", HLRZ, Jiilich, June 14-16, 1993. 



2. Introduction 



2.1. QCD — a challenge m computational physics 

Quantum chromodynamics (QCD) as the fundamental theory of strong interactions 
is closely related to quantum electrodynamics (QED), both being built on the prin- 
ciple of local gauge invariance. The gauge group of QCD is SU(3) and refers to an 
internal charge space of three "colour" degrees of freedom. 

Hadrons like the proton are composed of quarks bound by the exchange of 
gluons, the QCD analogs of the photons. The QCD action S thus has a very 
appealing (since geometrical) structure, involving gluon self-couplings which make 
it so different from QED. The proper evaluation of QCD, nevertheless, still presents 
a major challenge to the theorists: their goal still is to derive the spectrum and the 
matrix elements of hadrons from the basic QCD action. Such ab initio calculations 
would (a) imply a clean empirical verification of the theory, which could thus be 
tested against a huge set of empirical low energy phenomena and (b) allow to 
unfold QCD effects from experimental data in future studies of the standard model 
of electroweak interactions. 

Perturbation theory is obviously not adequate for the QCD evaluation within 
the low energy regime, because of the strong coupling involved; yet is has been 
proven to be applicable in the high energy limit of the theory, because of asymptotic 
freedom. 

In order to go beyond perturbation theory, we remember that quantum field 
theories can be, through their path integral formulation, cast into a statistical me- 
chanics problem. For this purpose, one has to continue from real to imaginary time, 
t —it. This Euclidean trick converts the phase factor exp(— iW) to a given "path" 
(field configuration) into a positive real weight, that can be readily interpreted as 
classical Boltzmann weight. So the technique is to estimate the (functional) path- 
integrals by generating a representative ensemble of field configurations $ with 
Boltzmann distribution 

P(<I>) ~ exp(-5[<I>]) (1) 

by a stochastic process and to "measure" observables within this ensemble. After 
space-time discretization this amounts to Monte-Carlo summation of high dimen- 
sional integrals (dimension at least 10®). Discretization implies a lattice regulariza- 
tion - and thus a proper definition - of the original continuum quantum field theory. 
The lattice formulation of QCD has been introduced by Kenneth Wilson in 1974^ 
and leads us directly into the arena of high performance scientific supercomputing^ 
(HPSC). Normally, one chooses a hypercubic four-dimensional lattice. This lends 
itself easily to parallel algorithms, at least when the updating within the Markov 
chain is carried out locally. For more details about our actual implementation on 
the Connection Machine CM-2, see Ref. 3. 



2.2. The Wupperial parallel computing setup 

Hopefully you can appreciate by now, that the challenge of "solving" QCD provides 
sufficient motivation to get involved in HPSC. Indeed, during the past years lat- 
tice gauge theory (LGT) has been an important motor for a variety of pioneering 
activities in parallel computing all over the world. In the University of Wuppertal 
the LGT group triggered four research groups from theoretical and experimental 
physics, electrical engineering and applied computer science into a joint initiative 
to launch an interdiscplinary computer laboratory for massively parallel computing 
in science and engineering. 

In October 1990, after two years of struggling with the funding agencies, the 
first Connection Machine in a German university (8K CM-2 with 256 MByte mem- 
ory, and a 10 GByte DATAVAULT funded by the Deutsche Forschungsgemeinschaft) 
was installed in Wuppertal and given into the hands of our students. Our opera- 
tional style is clearly 

Jit research-oriented Jit few-group Jit semi-dedicated Jit educational Jit 

The Connection Machines are run both by theoretical physicists and computer 
scientists in an informal (i.e. non computer center) style. We have developed a very 
nice cooperation with Thinking Machines Corporation in Boston, the manufacturer 
of our machines. We organize a joint seminar among the university groups involved 
and we are open for outside users to test their implementations on our machines. 

We would like to stress, that the availability of a "home" supercomputer is 
of tremendous importance to the research students in computational physics: they 
need a supercomputer to "touch" , on which they can try out their own ideas about 
algorithms, without an overly formal admission scheme common to supercomputer 
centers. Yet they have the compute power of a CRAY YMP at their disposal. In 
this manner, a university can fulfill its mission to challenge and train the creativity 
of a young generation of researchers going into computational science, which, after 
all, has been declared one of the goals of the US and EC grand challenge programs 
in HPSC. 

You recognize that we are highly commited to promote HPSC in the university 
environment, after its previous exodus into the supercomputer centers. This is 
a realistic strategy, due to the advent of low price, parallel supercomputers on 
the market. So today, we are facing a unique chance to boost the computational 
approach to scientific problems by the appropriate university training. 

In our case university and government agencies agreed that our pilot project 
should evolve with CM hardware. In this spirit we purchased a Connection Machine 
CM-5 (128 vector units, 1 GByte memory, 24 GByte scalable disk array), which at 
the time (October 1992) was the first of its kind in Germany. This pilot machine is 
presently performing 800 Mfiops in CMFortran on a typical QCD code like Wilson 
fermion propagators. We expect this performance to increase at least by a factor 
1.5 within the next half year. 



3. The Creutz experiment 

We think it is fair to say that Creutz' 1979 computer experiment ■* to "measure" 
the static qq potential opened a new era in computational hadron physics, based on 
the methods outlined in section . He produced first evidence that lattice techniques 
provide a sensible way to tackle nonperturbative effects in strong interaction physics. 
His conclusions in essence have been that quarks are confined, i.e. that the static 
potential exhibits a linear rise at large distances, and - just as important - that 
lattice results seem relevant to continuum physics. 

Let us quickly explain how to "measure" the potential on the lattice. We start 
from a closed rectangular path C(R, T) with extension R x T. The potential is 
extracted from so called Wilson loops W(R, T), which are defined to be the trace of 
path-ordered products of link variables Ufi(n) (associated to the link emanating in 
direction /j, from lattice site n) along such a path C(R,T). This loop construction 
corresponds to the world lines of a quark-antiquark pair at rest, separated by dis- 
tance R from each other and "travelling" over a time separation T. In Euclidean 
time, this observable will reveal the static "ground state" energy for asymptotic, 
i.e. for sufficiently large T values: 



So we analyse in the interval T > Tmm , with some reasonable cutoff Tmm and make 
use of the local mass 



as an estimator for the potential V^(R). For more details, we refer to Refs. 5,3,6. 

The crucial question in the context of lattice computations is: do we see contin- 
uum physics? Or put it differently: do we observe asymptotic scaling? The answer 
to this question is connected to the dependence of the gauge coupling parameter, 
g{a), on the lattice spacing, a, on the hypercubic lattice. In perturbation theory, 
this gauge coupling parameter can be shown to vanish in the limit a ^ 0, a feature 
which is called asymptotic freedom. Renormalization group theory yields in the 
two-loop approximation for the inverse function a{g): 



which form contains an energy scale, of the theory, such that the left hand side 
to this equation is dimensionless. The first two coefficients of the weak coupling ex- 
pansion of the SU (N) Callan-Symanzik /3-function are determined in perturbation 
theory to be 



VF(R, T) = 0(R) exp {-V{IL)T} . 



(2) 
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The number of colours in our case is # = 3. 
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Table. 1. The average plaquette action Sa, ^/a in lattice units, and the ratio ^/a/ Aj^, 
calculated by use of the perturbative approximation Eq. 4 from the bare /? and an effective 
coupling Pe = "^/Sa, respectively. 

Note that these formulas pertain to the so called zero-flavor (quenched) sector 
of the theory, which neglects fermion-antifermion loops, in the spirit of the valence 
approximation. In these relations, g stands for the so called bare lattice coupling, as 
it appears in the lattice action S. The lattice spacing a is determined by matching 
any dimensionful empirical quantity to its respective lattice result, at a chosen g. 
With a determined in this way, Eq. 4 can be used to estimate with different 
values of g. Asymptotic scaling is established, once the scale parameter A^ is 
verified not to depend on g any more. This is expected to occur at "sufficiently" 
small values of g. A priori, it is not clear, however, whether and where this regime 
can be reached in practical simulations, as we require, for the present purposes, a 
minimum physical lattice extent (about 1 fm) to avoid large finite size effects in the 
calculation, which amounts to an increasing number of lattice sites as a ^ 0. 

4. Physics discussion 

We are interested in the static qq potential, which is expected to be dominated by 
a Coulomb-like term at short distances with an R dependent coupling 

ViR)~Vo-aviaR)/R for i? small (6) 

and by a linearly rising term at large distances 

V(R) ~ KR for R large. (7) 

The coefficient K is related to the empirical string tension a by Ka~^ = a = 
(440MeT/)^. A " measurement" of K on the lattice thus provides the scale: the 
lattice spacing a(g) = \J u j K{g) to a given value of the bare lattice coupling g. On 
the other hand, one might be tempted to extract from the short-distance behaviour 
of the potential the strength parameter 



av{aK) = gl{aK)l{\^). 



(8) 
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Fig. 1. The ratio C = ^/a/ Aj^, in the two-loop weak coupling approximation Eq. 4, 
calculated by use of the bare coupling = (>/ jS and the /9_B-scheme g\ = 3Sn- The error 
band denotes the value C = 1.83 ± 0.07, extracted from the running of the interquark 
coupling aqq. 



The renormalized coupling gv can be connected to the bare lattice coupling g, and 
is due to run towards zero in a logarithmic fashion. In order to reduce the number 
of free parameters in this analysis, we investigate the force F(R) = —dV(R)/dR = 
—aqq/R"^ rather than the potential itself. The quantity aqq in the two-loop approx- 
imation has the form 

aqq{aR)= ^[bo\n{RaAR)-'' + bi/bo\n\n{RaAR)-''^ \ (9) 

where perturbation theory relates the scales Kr and A^. 

Physical masses like the quantity = can be retrieved from the lattice 
as inverse correlation lengths, i.e. their measurement asks for large spatial distances 
(in units of a) on top of the requirement of sufficiently large time separations T . So 
we must ascertain the inequality 

a<^<i.a (10) 

in order to attain continuum physics. In our simulations La was chosen to be larger 
than 1 fm. 

The string tension "measurement" should be performed at large values of R: 
K = —Yhhr^co F{R). Practically, this requires large lattice volumes and huge 
statistics. Note however, that in reality the extraction of a from the charmonium 



spectrum is based on a length scale of ~ 1 fm^! In both instances, one cannot 
neglect l/R corrections to the long distance potential, Eq. 7. For this reason, it 
is practically impossible to achieve complete parametrization independence of the 
string tension analysis, in the available range of qq separations. 




R 

Fig. 2. The reconstructed continuum potential V(R) (Eq. 15), measured at /? = 6.4 as a 
function of R with the corresponding fit curve. 

In the following we will present two different approaches to a potential analysis 
of our data. Version I is based on the Cornell parametrization commonly used in 
the charmonium analysis, while version II aims at extracting continuum information 
from short distance data as well, where lattice artifacts have not yet passed away. 

Version I. In the spirit of the potential analysis of quarkonia states, we start 
out from the Cornell ansatz for the potential^ 

V{R) = Vo + KR- e/R. (11) 

As described above, the fit range in R has to be restricted. Since we are aiming 
at the parameter K that governs the large i?-behaviour, we limit ourselves to the 
physical region aR > 0.3 fm, in which the experimental interquark potential has 
been reconstructed from the quarkonia states. We find that the last contribution 
to Eq. 11 can only be determined accurately for small (/-values, j3 = 6/(/^ > 6.0+. 

^Remember, that in the real world, with dynamical fermions, the confining part of the potential a 
la Eq. 7 will be screened by string breaking. 

+ For /? < 6.0, we have fixed e to the value determined from a fit at /? = 6.2, where we found 
e = 0.295 ± 0.017 
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Fig. 3. Comparison between deviations of the potential at /? = 6.4 from the continuum 
symmetry parametrization V(R) = Vo + K R - e/ R + f / , and 47rGi, (R) - (Eq. 13), 
which has been used to parametrize lattice artifacts. 

For a scaling analysis, it is very annoying to input data from different sources, 
as they are subject to different systematic errors. In order to keep control on the 
latter, we have carried out a comprehensive study using the analysis method just 
described over a broad range of /3: 5.5 < /? < 6.8. As a result, in Tab. 1 we 
can present the most accurate data set related to a dimensionful quantity ever 
determined in SU (3) lattice gauge theory. 

In Fig. 1 the combination C = -yo^/A-^jj-^r is plotted vs. /3, where the perturbative 
relation = 28.81Ai has been used. C has been computed via Eq. 4 (a) from 

the bare coupling, 28.81C = \/k'/ fig) (circles) and (b) from an effective coupling, 
13.88C = ^\/KI figs) (squares) as described in Ref. 5. The error band refers to 
the corresponding result from the running coupling analysis presented below. We 
conclude, that the strong scaling violations observed in the bare lattice coupling 
can be considerably reduced when passing to an effective coupling scheme. 

Version II. Here, we start from a parametrization® that accounts for lattice 
artifacts and thus enables us to make use of (more or less) the whole i?-region: 

T/(R) = T/o + A'i?-e fi^ + /47rGi(R)') + (12) 
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Our analysis method has been explained in more detail in a recent paper.® The 



lattice propagator for one-gluon exchange 

r d'k cos(kR) 

- (^4E,sm^(^./2) ^''^ 

occurs inside a correction term 

(5y(R) = e/(47rGL(R) - 1/i?) , (14) 
which is supposed to render the lattice-potential V(IV) rotationally invariant: 

V(R) = V(Il) + 6V(Il). (15) 
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Fig. 4. The on-axis interquark force at /? = 6.4 with different lattice artifact corrections. 
I = denotes the case without any correction, I = 1 holds for a pure lattice one-gluon- 
exchange correction, h denotes the "best" correction parameter obtained from the five 
parameter fit Eq. 12, and h is obtained from a constrained four parameter fit with / = 0. 
The horizontal line (with dashed error) is the asymptotic value —K for the force, calculated 
by use of the Cornell parametrization Eq. 11. 



This approach is remarkably successful as can be seen in Fig. 2, where the 
"corrected" potential values V(R) are plotted together with the corresponding fit 
curve. One might worry, how well the ansatz, Eq. 14, really covers the lattice 
artifacts. To answer this question, we compare in Fig. 3 the prediction (squares) of 
Eq. 14 against 6V(lV)/(el) (circles), as computed from the data. We find, that for 
R > the correction procedure works very well.§ 

^The fit incorporates five free parameters on a data set of 70 points. Tlie fit range lias been 
restricted to i? > \/3. 



We consider next the computation of the running coupling aqq(aR). We start 
out from the corrected potential data and construct the force, in order to eliminate 
the Vo contribution. The numerical differentiation is carried out in such a way that 
a pure Coulomb potential would yield a constant coupling. This is a reasonable 
approach to the physical situation, where the coupling is expected to vary only 
weakly with R. 

We emphasize that during this stage, i.e. in the computation of the force from 
the corrected potential, we are independent of the previous parametrization Eq. 12, 
as we perform the numerical differentiation (pointwise) on the very potential data. 
The parametrization has only been utilized to determine the lattice corrections. 




Fig. 5. The running coupling aqq in units of the string tension together with a one 
parameter fit with the two-loop formula, Eq. 9. 

In Fig. 4 we display the force for the on-axis separations, where the difference 
between lattice (Eq. 13) and continuum propagators is most pronounced. The plot 
is meant to illustrate the rather weak sensitivity of the force data with respect to 
a reasonable change in the quantity el controlling our lattice correction. We show 
the cases / = (pure Coulomb, triangles), / = 1 (pure lattice one-gluon exchange, 
diamonds), and li = 0.639(35) = 0.216(ll)/ei (circles), the latter being our optimal 
value. One might worry about the impact of the term on the lattice-corrected 

force. For this reason, we also tried a fit with / constrained to the value / = 0; 
as a result we find the value I2 = 0.600(58) = 0.170(16)/e2 (squares). The figure 
moreover contains the asymptotic value of the force, —K, as obtained from the 
Cornell approach (horizontal error band). 

The difference between the li and I2 results will be interpreted as systematic 



uncertainty and is included into the error of aqq(aR). We note, that the R = 2 
value of the force has not been used in the further analysis, as it involves using 
the potential data on i? = 1 or \/2, which might still be polluted by discretization 
errors. 

We show the results® of the aqq analysis in Fig. 5. The curve refers to a fit 
to the two-loop formula Eq. 9, with Ar = 0.572(22)-^^. It can be seen that the 
short distance data on our lattices reveal the running of the coupling expected from 
perturbation theory. A comparison of the scale Ajjg- calculated from Ar and from 
Ai is contained in Fig. 1 and shows beautiful consistency of the short and long 
distance features. 

5. Conclusion 

It appears that we are finally arriving at a point where, with the available lattice 
techniques, we can make contact to the perturbative weak coupling regime, i.e. to 
continuum physics - at least in the quenched sector to QCD. This demonstrates 
indeed that parallel supercomputers like the Connection Machine are useful research 
tools in HPSC! 
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